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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 
AN ANALYTICAL STUDY OF HEAT REQUIREMENTS FOR ICING 
PROTECTION OF RADOMES 


By James Р. Lewis 


SUMMARY 


The heat requirements for the icing protection of two radome con- 
figurations have been studied over a range of design icing conditions. 
Both the protection limits of a typical thermal protection system and 
the relative effects of the various icing variables have been determined. 
For full evaporation of all impinging water, an effective heat density 
of 14 watts per square inch wes required. When a combination of the 
full evaporation and running-wet surface systems was employed, a heat 
requirement of 5 watts per square inch provided protection at severe 
icing and operating conditions. 


INTRODUCTION 


Radar is becoming of increasing importance in the design and 
operation of aircraft. Successful operation of these aircraft demands 
that the performance of the radar system including.the radome be unim- 
paired by environmental factors including icing conditions. Recent 
experimental investigations have shown that radomes mounted in the nose 
of an aircraft are very susceptible to icing and that this radome icing 
has serious effects on the radar operation resulting from a marked 
decrease in the transmission efficiency and a deflection of the radar 
beam. Protection of the radome against icing, therefore, is required. 
This protection can be achieved in several ways, including: applying 
a temperature-depressant material to the radome surface, or heating the 
radome surface sufficiently to prevent the formation of ice. This report 
will consider only the thermal protection method. 


The determination of the heating requirements and the performance 
of a thermal icing protection system requires many complex calculations. 
Much of the basic information required for such а calculation, particu- 
larly for a body such as а radome operating at high speeds, is not always 
readily available nor in a form directly applicable for engineering 
design purposes. For these reasons, the performance of a radome thermal 
icing protection system has been studied at assumed operational and icing 
conditions. The objectives of this study, which was conducted at the . 
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МАСА Lewis laboratory, were to determine the protection requirements 
for a typical thermal icing protection system and to study the relative 
effects of the various icing variables on these protection requirements. 


ANALYSIS 


Thig investigation is divided into two parts: (1) а study of pro- 
tection requirements for complete evaporation of all the water impinging 
upon the radome, resulting in а dry surface, and (2) a study of protec- 
tion requirements for the case of & running-wet condition, that is, only 
sufficient heat is supplied to maintain the coldest point on the surface 
of the radome at 32° F and thus prevent the formation of ice on the 
radome. Іп this latter method, the excess water not evaporated from 
the radome surface will flow aft and freeze on the unheated portions of 
the aircraft. 


Evaporation of Impinging Water 


For the case of complete evaporation of all the impinging water, 
the heat- and mass-transfer relations given. in references 1 and 2 may 
be written in the following form i 


Hay(tg-to) + 2.82 ІҺ,/К(е,/рі - ео/ро) + Cyilgy(tg-to) 
(1) 
PE os 


All symbols are defined in appendix А. In order to simplify the caleu- 
lations, average values are used in this equation rather than calcu- 
lations made at local points all over the radome surface and integra- 
tion of these results. In addition, the heat quantities are the heat 
requirements at the radome surface and do not represent the heat-source 
requirement. 


The heat required for evaporation is 


de = Mayl (2) 


The heat available for evaporation is . 


de = 2.82 hy L(e,/p1 - ео/ро) (3) 
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For complete evaporation equation (2) must equal equation (3), or 


2.82 (е„/ру - ер/Ро) = шау/Вау (4) 
The average rate of water impingement is given by 
m y = 0.3296 Ел Vo Ap/Ag w (5) 


The collection efficiency E, is defined as the ratio of the 
amount of water actually impinging upon a body to the amount of cloud 
water which would be swept out by the area of the body projected in the 
Plight direction. This collection efficiency is a function of body 
size and shape, airspeed, temperature, pressure, and water droplet size. 
The values of collection efficiency used in this analysis were obtained 
from the data of reference 3 and unpublished experimental data. 


The latent heat of vaporization L was taken as 1060 Btu per 
pound, К as 1.0, V ав 0.87 Vg, and г ав 0.85. Assuming values 
of % and corresponding values of po and Ең for a given condition 
of airspeed and effective heat input а, equations (1), (4), and (5) 
were solved simultaneously by trial and error for the limiting values of 
liquid-weter content. This value of liquid-water content represents 
the maximum value for the assumed conditions for which all the impinging 
water will be evaporated within the heated surface area of the radome. 
An airspeed of 600 miles per hour was assumed throughout the study 
together with an effective heat density of 2100 Btu per hour per square 


foot, equal to 47 watts per square inch. This heat density was con- 


sidered to be uniform over the surface area of the radome. 


Running-Wet Condition 


For the case in which the coldest point on the radome is just 
529 Р апа no ice forms ; the following solution was used: 


The local rate of impingement is given by 


m = 0.3296 VgwB (6) 


where the local collection efficiency В was determined from the data 
of reference 3 and unpublished experimental data. 
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By use of the relative heat factor b = сұ m/h, equation (1) for 
this condition may be written 


(32-to) (145) + 2:82 IK(0.18 - ео) 


ро А21/20 


(7) 


For this case, L was taken as 1075 Btu per pound. Preliminary 
calculations at various points along the radome surface indicated that 
the coldest region on the radome occurs at the rearmost point on the 
radome. The limit of impingement was also found to be very close to 
this rearmost point. This results in a value of the wetness factor K 
of 1.0. As in the case for complete evaporation, equations (6) and (7) 
were solved for the limiting value of liquid-water content for various 
air temperatures at the assumed airspeed and effective heat density 
values. In addition, the required value of liquid-water content for a 
running-wet surface for the case of an unheated radome was obtained by 
setting the q/h term in equation (7) equal to zero. 


Icing Conditions 


The degree of protection afforded a vulnerable aircraft component 
by a thermal icing protection system is dependent upon the icing condi- 
tions that will be encountered as well as upon the availability of heat 
and the system efficiency. Thus, in the design and in the appraisal of 
a protection system, the expected icing conditions must be studied and. 
established in order to obtain answers that are of reasonable engineer- 
ing validity. The important variables that must be considered in 
defining an icing condition are the cloud liquid-water content, the 
water droplet size and size distribution, the air temperature and pres- 
sure, and the extent and frequency of occurrence of a particular type 
of cloud. Extensive studies of these factors and their combinations 
have been made by the NACA. Statistical studies of icing conditions 
and methods of determining the proper combination of the important icing 
variables have been reported in references 4 and 5. 


Tn the selection of the variables defining an icing condition, the 
collection efficiency of the body must also be considered since it is 
as important as the cloud water content and cloud extent in determining 
the severity of a particular icing condition. The collection efficiency 
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of the radomes considered in this analysis was used (as indicated in 
ref. 4) to determine the particular combination of values of liquid- 
water content and droplet size for a particular frequency of occurrence 
that resulted in the maximum rate of impingement. 


Since the heat requirements for a thermal protection system are 
dependent on the air temperature, it was decided to employ a temperature- 
altitude relation representative of icing conditions. Reference 4 pre- 
sents the observed variation of air temperature and pressure altitude 
in icing conditions. Considerable variation in air temperature for a 
given altitude is shown by these data; and for this reason and also to 
obtain a more realistic basis for appraisal of the protection system, 
the temperature-altitude curves shown in figure 1 were selected. One 
curve represents the average of the data of reference 4, while the 
second curve is a fairing through the points of lowest temperature 
reported in reference 4. The NACA standard atmosphere is also given 
in figure 1 for comparative purposes. From these temperature-altitude’ 
relations and from an assumed exceedence probability of 1 in 1000, the 
particular combinations of liguid-water content and droplet size which 
gave maximum rate of impingement were chosen from the curves of refer- 
ence 5. The liguid-water content and droplet size corresponding to 
cumulus clouds were taken for the low-temperature condition ав repre- 
senting extremely severe icing conditions, while layer cloud values 
were taken for the average temperature conditions as being typical of 
average icing conditions, especially with respect to extent of the con- 
dition апа ав the limit for which the full evaporation system would 
provide protection. The resultant curves of liquid-water content and 
droplet size against air temperature are shown in figure 2. The drop- 
let size distribution assumed in the analysis is shown in figure 3. 


Miscellaneous Ássumptions 


in addition to the icing conditions, several other factors were 
assumed for the purposes of the analysis. Two radome configurations 
were investigated, half-sections of which аге shown in figure 4, 
together with the pertinent dimensions. Both radomes, which were 
assumed to be nose installations, were portions of ellipsoidal bodies 
of revolution. It was assumed that protection was required to the rear 
of the radomes. | 


The convective heat-transfer coefficient was calculated from 
unpublished experimental data obtained from tests of similar bodies 
in the icing research tunnel. 


6 NACA RM Е55А22 
RESULTS AND DISCUSSION 
The results of this analysis are presented for the two radome con- si 


figurations for the assumed icing conditions. 


A-Radome 


1422 


Low-temperature cumulus condition. - The performance of the pro- 
tection system for the blunt A-radome at the low-temperature  cumulus- 
cloud condition is shown in figure 5. The assumed icing conditions | 
taken from figure 2 are designated as the icing limit, and it is assumed 
that no protection is necessary for conditions to the left of this curve. 
For the case in which no heat is applied te the radome, an ice-free 
running-wet surface for certain conditions results from the kinetic tem- 
perature rise. This condition is obtained for all temperatures above 
approximately 2° Р (equivalent to approximately 8000 ft). With the 


assumed heat density of 42 watts per square inch, the limit for the 


ice-free running-wet condition becomes approximately -18° F. For all « 
air temperatures to the right of this curve and for the values of 
liquid-water content shown, the surface temperature will be equal to 

or greater than 32° F with varying percentages of the impinging water 
being evaporated. Ав indicated by both the heated and unheated running- 
wet curves, the requirement for this condition is almost independent of 
variations in liquid-water content and is almost entirely dependent upon 
the air temperature. 


For the case of full evaporation, in which the radome surface is 
maintained dry, the limiting liquid-water content varies from approxi- 
mately 0.45 gram per cubic meter at -30° F to 0.5 gram per cubic meter 
at 209 F. The surface wlll be kept dry for all liquid-water and air- 
temperature conditions below the full evaporation curve. For the case 
of full evaporation, the calculated average surface temperature varied 
from approximately 659 to 769 F. In contrast to the case of the running- 
wet condition, the requirement for full evaporation is almost independent 
of variations in air temperature. 


The area below the icing-condition-limit curve that is not protected 
either by full evaporation or the running-wet condition is seen to be 
rather small. Іп this region, ice resulting from both direct water 
impingement and from runback and refreezing on the radome surface will 
be obtained. 


In order to obtain an estimate of the heat density required for 
protection over the full range of expected icing conditions, calculations | ^ 
were made of the variation of the heat requirement with water content 
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and with air temperature for the full evaporation апа running-wet sys- 
tems, respectively, at specific values of altitude, droplet size, tem- 
perature, and liquid-weter content. These results are presented in 
figures 6 and 7 for the A-radome. The effect of kinetic heating is 
indicated in figure 6 by the evaporation of water up to 0.09 gram per 
cubic meter without the application of heat. The heat requirement for 
full evaporation is seen to vary almost linearly with liquid-water con- 
tent; the heat requirement for а running-wet surface also approaches 

а linear relation with air temperature. From an extrapolation of the 
relations of figure 6, it is determined that protection by means of full 
evaporation alone over the full range of icing conditions would require 
an effective heat density of approximately 14 watts per square inch. 
Full protection by a combination of the evaporation and running-wet 
surface systems over the entire range of expected. icing conditions could 
be achieved with an effective heat density of approximately 5 watts per 
square inch. 


Average—temperature layer-cloud condition. - The performance of 
the thermal protection system for the A-radome at the average- 
temperature layer-cloud condition is shown in figure 8. In this case 


the assumed effective heat density of 47 watts per square inch is suf- 


ficient to provide protection over the entire range of icing conditions. 
Full evaporation of the impinging water is obtained up to approximately 
0.45 gram per cubic meter corresponding to an altitude of approximately 
9800 feet and a temperature of 12° F. A heat density of 6.5 watts per 
square inch would provide full evaporation over the whole range of icing 
conditions. The minimum heat density that would provide protection by 
a combination of systems (full evaporation and running-wet surface) is 
approximately 2.5 watts per square inch. The limit for the no-heat 
running-wet surface condition is 5.59 Р as compared with 2° F for the 
low-temperature cumulus-cloud condition (fig. 5). 


B-Radome 


Similar results for the narrow B-radome are presented in figure 9 
for the low-temperature condition. The limits of protection for the 
running-wet surface condition both with and without heat are approxi- 
mately 2° F less than for the blunt A-radome. The average limiting 
liquid-water content for full evaporation is approximately 0.53 gram 
per cubic meter as against 0.48 gram per cubic meter for the A-radome. 
Thus, despite a considerable difference in both fineness ratio and size 
of the two radome configurations, the performance of the protection 
systems is almost the same. For this reason the results for the 
B-radome at the average-temperature layer-cloud condition have been 
omitted. 
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CONCLUDING REMARKS 


In the analysis presented herein it was necessary to make several 
assumptions. It is believed that most of these assumptions, including 
the icing conditions, are of reasonable validity. The most important 
assumptions that might be questioned are the impingement efficiency and 
the assumption of an average or uniform effective heat density. The 
impingement efficiencies used. in the analysis are based upon the values 
for spheres given in reference 3 and upon experimental results for віш- 
ilar radomes obtained at lower airspeeds. For the case of the running- 
wet surface, the value of the assumed local impingement efficiency is 
relatively unimportant since, as shown by the results of the analysis, 
the heat requirement is almost independent of the amount of water caught. 
For the case of full evaporation the heat requirement, as indicated by 
the results of figure 6, is directly dependent on the amount of water 
caught. The total collection efficiencies are regarded as accurate 
Within at least +10 percent. Based upon the results of figure 6, the 
heat requirements in the range of interest have approximately the same 
degree of accuracy. 


А more important limitation of the results for the case of full 
evaporation is the fact that a uniform heat density over the radome 
surface was assumed and the calculations were made on an average basis 
rather than by computing the system performance on a point-to-point 
basis from the radome nose aft. The attainment of both a uniform heat 
density and surface temperature is impossible and, practically, even 
the attainment of uniform heat density would be extremely difficult. 
It 1s believed, however, that the use of average values in the calcu- 
lation of the performance of the protection system is valid for the 
purposes of this analysis as indicating within the limits of engineering 
accuracy the order of magnitude of the limits of performance of the 
protection system and the variation of these limits with the important 
icing variables. 


The results of the analysis have indicated that icing protection 
of a radome by a thermal protection system can be achieved with reason- 
able values of heat density even at extreme combinations of operating 
and icing conditions. For full evaporation of all impinging water, an 
effective heat density of 14 watts per square inch is required. By 
employing a running-wet surface system over part of the temperature 
range and full evaporation over the remainder, protection over the full 
range of icing conditions can be achieved with an effective heat density 
of 5 watts per square inch. The heat requirement for full evaporation 
is dependent primarily upon the rate of water catch or, in terms of the 
icing condition, the cloud liguid-water content. For the running-wet 
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surface condition, the heat requirement is primarily a function of the 
ambient air temperature. The use of the running-wet surface system will 
be dependent not only on the effects of a water film on radar operation 
but also on the tolerance of the aircraft for runback ice formations aft 
of the radome. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX - SYMBOLS 


The following symbols are used in this report: 


А area, sq ft 

b relative heat factor, me /h, dimensionless 

с, specific heat of air at constant pressure, 0.24 Btu/(1b) (°F) 
es specific heat of liquid water, 1.0 Btu/(1b) (°F) 

а. droplet diameter, microns 

En collection efficiency, dimensionless 

e partial pressure of water vapor, (corresponding to saturated 


air), in. Hg 


g acceleration due to gravity, 32.2 ft/sec? 

h dry convective heat-transfer coefficient, Btu/(hr)(sq ft) (°F) 
J mechanical equivalent of heat, 778 (ft)(1b)/Btu 

K surface wetness fraction, dimensionless 

L latent heat of vaporization of water, Btu/lb 

m rate of interception of water, lb/(hr)(sq ft) 

p absolute static pressure, in. Hg 

а effective surface heat density, Btu/(hr)(sq ft) or watts/sq in. 
r kinetic energy recovery factor, dimensionless 

t temperature, ОҒ 

ү air velocity, ft/sec 

W cloud liquid-water content, g/cu m 


В local collection efficiency, dimensionless 
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Subscripts: 

ау average value 

e evaporation value 

p projected. in stream direction 
S surface value 

O stream value 

1 local value 
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Figure 1. - Assumed temperature-pressure relations for icing 
conditions. 
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Figure 2. - Variation of cloud liquid-water content and droplet size with 
air temperature. 
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Figure 5. 
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Figure 4. - Assumed radome configurations. 
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Figure 6. - Variation of heat required for full evaporation with liquid- 
water content. A-Radome; airspeed, 600 miles per hour; pressure altitude, 
15,000 feet; temperature, -11.5° F; droplet diameter, 20 microns. 
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Figure 7. - Variation of heat required for running-wet 
surface with air temperature.  A-Radome; airspeed, 600 miles 
per hour; pressure altitude, 20,000 feet. 
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Figure 8. - Performance of protection system at average-temperature » layer-cloud conditions. 
A-Radome; airspeed, 600 miles per hour; effective power density, 2100 Btu per hour per 


square foot, or 4 watts per square inch. 
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Figure 9. - Performance of protection sys at low-temperature, cumulus- 
cloud condition.  B-Radome; airspeed, 600 miles per hour; effective power 
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density, 2100 Btu per hour per square foot, or 4 Watts per square inch. 


